Abstract. The interpretation of high proper motion 
white dwarfs detected by Oppenheimer et al. (2001) was 
the start of a hvely controversy. While the discoverers 
identify a large fraction of their findings as dark halo mem- 
bers, others interpret the same sample as essentially made 
of disc and/or thick disc stars. We use the comprehen- 
sive description of Galactic stellar populations provided 
by the "Besangon" model to produce a realistic simulation 
of Oppenheimer et al. data, including all observational se- 
lections and calibration biases. The conclusion is unam- 
biguous: Thick disc white dwarfs resulting from ordinary 
' hypotheses on the local density and kinematics are suffi- 
p\| . cient to explain the observed objects, there is no need for 
Q^" halo white dwarfs. This conclusion is robust to reasonable 
, changes in model ingredients. The main cause of the mis- 
^0 ' interpretation seems to be that the velocity distribution of 
i: a proper motion selected star sample is severely biased in 
,—1 favour of high velocities. This has been neglected in previ- 

■ ous analyses. Obviously this does not prove that no such 

■ objects like halo white dwarfs can exist, but Oppenheimer 
et al. observations drive their possible contribution in the 
dark matter halo down to an extremely low fraction. 
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1. Introduction 



In a recent paper, Oppenheimer et al. (2001 ) (hereafter 
OHDHS) give evidence for a large number of faint white 
dwarfs detected in digitized photographic plates from the 
SuperCOSMOS Sky Survey^. Their interpretation of these 
stars as dark halo members was the start of a contro- 
versy. While iRcid et al. (200l| ) and more recently |Grafi] 



are a function of age, and constrained by the Galactic 
potential through the Boltzmann equation. The average 
kinematics of white dwarfs in a volume limited sample 
in the solar neighbourhood result in a velocity ellipsoid 
(42.1,27.2,17.2) and an asymmetric drift of 16.6 km s^^. 

The thick disc population is modeled as originating 
from a single epoch of star formation. The adopted thick 



(2001 ) interpret OHDHS stars a s the high velocity tail of 



the disc/thick disc populations, Koopmans fc Blandford 



(2001 ) argue that there is a statistically significant break 



in the velocity distribution, which they interpret as the 
thick disc/halo break. They propose a dynamical mecha- 
nism to produce these "halo" objects out of the disc. An- 



other line of arguments has been developed by Gibson & 



Flynn (2001) who point out a number of biases in OHDHS 
analysis all adding up to produce an overestimation of the 
local density of halo white dwarfs, but the authors do not 
discuss their nature Tn addition, Hansen (2001 ) finds a 
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thin disc white dwarf population. In order to clarify the 
situation, we use the Besangon model of stellar popula- 
tions in the Milky Way to simulate OHDHS observations. 
In the next section we describe model hypothesis and the 
construction of the simulated sample, the third section 
being dedicated to the discussion and conclusions. 

2. Simulations 

We have performed simulations from the Besangon Galac- 
tic model using the following assessment: disc stars are 
assumed to have ages from to 10 Gyr, with appro- 



priate kinematics as a function of age (Haywood et al 
1997 ). The density laws are Einasto discs, w hich are close 
to sech^ vertically ( Bienayme et al., 1987a b|). The white 
dwarf luminosity function is from Sion fc Liebert (1977 ) 
and the photometry is derived from Bergeron et al. (1995 ) 
model atmospheres. The ellipticities of these components 
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disc density law is described in Elcylc & Robin (2001), 
as fitted to large sets of available star counts. The local 
density is 10"^ pc'^ or 7.1x10"'' Mopc"^ {My <8), that 
is about 6% of the disc local density. In order to ensure 
continuity of the density gradient at z = 0, the density law 
is made quadratic near the plane (z < 337 pc). Further 
out from the plane, the density law follows an exponential 
of scale height 850 pc. 

The white dwarf luminosity function (hereafter 
WDLF) of the thick disc has been computed by Garcia- 
Berro (private communication) from the model of [sern et 



al. (1998) assuming a Salpeter initial mass function (IMF) 



and an age of 12 Gyr. We attempted to normalise the 
WDLF by computing the number of white dwarfs rela- 
tive to main sequence stars, assuming that all stars with 
a mass greater than the mass at the turnoff (Mto = 0.83 
Mq) are now white dwarfs. We consider a two slope IMF, 
the lower mass parts being constrained by available star 



counts (Reyle fc Robin, 2001) accounting for a binary cor- 
rection: 



dN/dm oc m'^-^^; Mto <m < 8Mq 
dN/dm cx m"^ '^^; O.IMq < m < Mto 
The ratio of white dwarfs to main sequence stars is 
therefore: 



WD. 



0.74 " ■ 



SP 



= 0.43 



r TO m-o.75rf„j 

The factor 0.74 ensures the continuity at the turnoff. 
However, it appears that the predicted number of thick 
disc white dwarfs, with these assumptions, is much higher 
than the number of observed white dwarfs in the OHDHS 
sample. Several reasons could be invoked: First, we do 
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not know how the OHDHS sample is complete, specially 
near the limiting magnitude. Second, hypotheses which 
have led to the estimation of the thick disc WDLF have 
not been tested yet because these stars are not identi- 
fied as such in current white dwarfs samples: Uncertain- 
ties remain on the IMF slope (at high mass, with no ob- 
servational constraints, we have chosen the conservative 
Salpeter IMF, while at low mass an IMF slope a = 1, still 
within our error bars, drives down the fraction to 0.33). 
Uncertainties remain also on the star formation history, 
which may be more complex than a single burst. These 
reasons led us to normalize the thick disc WDLF in order 
to be in agreement with the OHDHS sample: the frac- 
tion of white dwarfs over sub-turnoff stars drops down to 
20%, about the value of Hansen (2001 ), corresponding to 



a local density of 5xl0~^ stars pc~^. This normalisation 
is also in agreement with the number of white dwarfs in 
the expanded LHS (Luyten Half Second proper motion 
catalogue) white dwarf sample (Liebert et al., 1999). The 
velocity ellipsoid of the thick disc is constrained by photo- 



astrometric survey in several Galactic directions (Ojha 
et al., .996) to be (67,51,42) while the asymmetric drift is 
53 km s^^. 

Dark halo white dwarfs are modeled assuming that the 
dark halo is filled with ancient white dwarfs with a certain 
factor /. The dark halo local d ensity is 8 x 10~'^ Mq pc" '^ 
from dynamical constraints ( Bienayme et al., 1987a| ). 
Hence the local density of halo white dwarfs is assumed 
to be / times this value. As a first guess compatible with 
microlensing constraints, we set / = 0.1 in this simula- 
tion. We use Bergeron et al. (1995) model atmospheres 
for hydrogen white dwarfs of mass 0.6 MqIo the turnover 
(My < 17), and [Chabricr et al. (2000[ ) models for the 
cooler part of the sequence. The halo white dwarfs lumi- 
nosity function has been computed by p^sern et al. (1998 ) 
with the ad hoc Chabricr (1999) IMF and an age of 12 or 
14 Gyr. The space velocities of the dark halo are computed 
assuming the same kinematics as for the spheroid: A ve- 
locity ellipsoid (131,106,85) and a null rotation, the LSR 
rotational velocity being 229 km s^^ (Bienayme et al., 
1987a|)l 

The colours are computed by transforming the stan- 
dard Johnson-Gousins system to the observed one using 
the equations given in the SuperCOSMOS Sky Survey de- 
scription: 

Bj = B -0.28 X (B-T/);for-0.1 <B-V <1.6 
R59F = i? - 0.006 - 0.059 x (i? - /) + 0.112 x (i? - 1)^ 
-f 0.0238 X (R-If 

The selection of the simulated sample has been done 
strictly following the OHDHS process. We simulate a field 
of 4165 deg^ at the south Galactic cap. We select a sam- 
ple of stars with R59F < 19.45, and proper motions in 
the range 0.33 to 3" yr~^, assuming a mean error of 
0.04" yr^^. In order to reproduce the space velocities de- 
duced by OHDHS, from the resulting sample, we compute 
the reduced proper motions, absolute magnitudes Mbj, 



distances and space velocities U and V using their formu- 
lae. 



3. Results and discussion 

Fig.|l|a reproduces the observed distribution of the 
OHDHS sample in the (U,V) plane. The ellipses indicate 
the la and 2a velocity dispersions as defined by OHDHS 
for the old disc (right) and the halo (left). Objects out of 
the 2a old disc ellipse are interpreted as dark halo white 
dwarfs by OHDHS (squares with boxes). This figure can 
be compared with the simulated distribution of the se- 
lected white dwarfs represented in fig.|l|b. The disc is coded 
with crosses, the thick disc with circles, and the dark halo 
with asterisks. Ellipses are the same as in figj^a for com- 
parison. 

The vast majority of OHDHS white dwarfs out of the 
2a old disc ellipse can in fact be interpreted as thick disc 
white dwarfs, as expected from ordinary hypotheses on 
the local density and kinematics of the thick disc white 
dwarfs, although it is not possible to be su re that none o f 
these white dwarfs belong to the dark halo. Hansen (2001 ) 
argues that the OHDHS white dwarfs are not old enough 
to account for the thick disc population, given that most 
of the white dwarfs in the thick disc should be fainter in 
the case of a single burst of star formation in the thick disc 
12 Gyr ago. However, the simulation shows that, even in 
the burst case, the number of younger white dwarfs in the 
thick disc remains large enough to explain the OHDHS 
sample. They come from stars for which the time spent 
on the main sequence is not negligible. Actually, very few 
of the simulated stars are beyond the turnover in the white 
dwarf colour-magnitude diagram. 

Another difference that appears in fig.|^ is the number 
of high velocity stars: Gonsidering an age of 12 Gyr for the 
halo white dwarfs and a fraction / = 0.1, 10 objects are 
expected to have V < —400 km s^^ (9 simulations were 
performed to compute a mean value for this number of 
high velocity objects). None of the OHDHS white dwarfs 
are present in this part of the diagram, indicating that 
the halo white dwarf local density is no more than 1% of 
the dark halo local density, in agreement with the values 



derived by Koopmans & Blandford (2001) and Gibson & 
Flynn (2001| ). However, this upper-limit remains uncer- 
tain as the completeness of the OHDHS sample is not 
precisely determined. Furthermore, if considering an halo 
age of 14 Gyr, the expected number is 2. 

The photometric estimation of the distances could in- 
troduce a bias. They are determined from a linear colour- 
magnitude relation, with uncertainties around 20%. The 
linear relation is no longer valid for old white dwarfs be- 
yond the turnover in the colour-magnitude diagram. Using 
this relation leads to an overestimation of their luminosi- 
ties and consequently, their distances and V velocities. 
However, very few stars on the OHDHS sample are af- 
fected by this bias (see figure 1 of Koopmans fc Blandforcj 



4 



Celine Reyle et al.: On high proper motion white dwarfs 



1 




-600 -400 -200 -600 -400 -200 

V (km/s) V (km/s) 



Fig. 1. Velocity space distribution of white dwarfs in the Oppenheimer et al. (2001) sample: a. Original observations. 
Objects with boxes are interpreted as halo stars, b. Simulation with a white dwarf local density being 10% (/ = 0.1) 
of the dark halo local density, and an age of 12 Gyr. Crosses: disc stars, circles: thick disc stars, asterisks: dark halo 
stars. 



(2001)). On the contrary, the simulated sample consider- 
ing a 12 Gyr dark halo with / = 0.1 contains some of these 
old white dwarfs, which are the stars one is looking for to 
explain the dark matter halo. 

Another strong bias comes from the selection of high 
proper motion stars which isolates a high velocity tail of 
each population. The ellipses in fig.|l| characterize kine- 
matics of a complete halo/thick disc sample, but are in- 
correct for a high proper motion sample. The simulated 
sample can be separated in three distinct populations and 
it well reproduces the observed dispersions and rotation 
lag. Stars within the la old disc ellipse are mainly disc 
stars, in contradiction with Koopmans & Blandford (2001) 
who neglected this population. Out of the 2a old disc el- 
lipse, stars are dominated by thick disc stars and a few 
disc stars. The fraction / = 0.1 for an halo of 12 Gyr 
clearly overestimates the number of observed stars in the 



Table 1. Rotational lag and velocity dispersion au for 
disc, thick disc and halo white dwarfs. Values in a complete 
sample are those in the original model. Values from the 
selected sample come from the simulation with the high 
proper motion and apparent magnitude selections. In the 
selected sample the computed velocity dispersion and lag 
are about twice what it would be in a complete sample. 







complete sample 


selected sample 




disc 


16.6 km s"' 


30 km s"' 


lag 


thick disc 


53 km 


122 km s"^ 




halo 


229 km s"^ 


336 km s~^ 




disc 


42.1 km 


64 km s"^ 




thick disc 


67 km s^^ 


112 km s"^ 




halo 


131 km s"^ 


207 km s"^ 



la halo ellipse and even more outside. 



larger asymmetric drift, such as 89 km s (as in Robin & 



As shown in table 1, due to the observational biases 
in the sample, the mean asymmetric drift of the resulting 
sample is much larger that the normal value for a complete 
sample in the solar neighbourhood. The expected values 
for halo white dwarfs are also given in table 1. There are 
much larger than the values found by OHDHS, implying 
that the halo population cannot dominate the sample. 

It should also be noted that, even if a large uncer- 
tainty remains on the true local density of thick disc white 
dwarfs, our conclusions are robust to changes of hypoth- 
esis within this uncertainty, since they are mainly based 
on the velocity characteristics of the sample. Uncertain- 
ties on the thick disc kinematics also exist. Adopting a 



Chen (1993)), would slightly displace the thick disc white 
dwarfs towards negative V, but our conclusions remain 
unchanged. Had we adopted an asymmetric drift small 
enough to displace most of the thick disc stars inside the 
2cr old disc ellipse, the fraction /=0.1 would underestimate 
the number of stars with —400 km s^^ <V< —200 km s~^, 
but a higher value for / would dramatically overestimates 
the number of high velocity stars with V< —400 km s~^. 

4. Conclusions 

The simulated sample, taking into account all observa- 
tional bias, gives an unambiguous result: Most high veloc- 
ity white dwarfs in the OHDHS sample can be safely in- 
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terpreted as disc and thick disc stars. The expected contri- 
bution of spheroid white dwarfs, assuming their mass den- 
sity to be 1.3x10^5 Mqpc-^ ( [Gould et al., 199S| ) would 
be 0.5 stars. The OHDHS sample has been misinterpreted 
because the bias introduced by the selection of high proper 
motions has been neglected. Although some of the stars 
in the sample may be part of the halo, it is not necessary 
to call for exotic objects such as white dwarfs in the dark 
halo. However, this sample provides a direct observation 
of thick disc white dwarfs. Such observations will lead to 
the construction of the luminosity function for thick disc 
white dwarfs, and hence to a better understanding of the 
Milky Way's history. 

The simulation also shows that a sample dominated 
by the dark halo white dwarfs would have much larger 
velocity dispersions than the OHDHS sample. The num- 
ber of expected candidates with a V velocity less than 
—400 km s~^, compared with the null detection, implies 
that the fraction of dark matter halo made of such objects 
is less than 1% for a 12 Gyr halo, or that the halo is older 
than 12 Gyr. 

In a more general way, this study shows that care 
must be taken when one tries to interpret the nature of 
high proper motion wh it e dwarfs in photog r aphic surveys 
(e.g. |Knox et al. (1999|) ; |lbata et al. (2000[ ); [Monet et af 
(200Cl(]J, as the thick disc population is a non-negligible 
part of high proper motion selected sample, and as thick 
disc white dwarfs cannot be easily distinguished from halo 
white dwarfs. As also shown by Creze et al. (2001 ), the 
search for halo white dwarfs needs to be performed in 
deeper astrometric surveys, when the thick disc becomes 
non dominant over the expected halo white dwarfs popu- 
lation in high proper motion selected samples. 
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